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Abstract 

We present a detailed analysis of the phase transition in the stan- 
dard model at finite temperature. Using an improved perturbation 
theory, where plasma masses are determined from a set of one-loop 
gap equations, we evaluate the effective potential Veff{^p,T) in next- 
to-leading order, i.e., including terms cubic in the gauge coupling g, 
the scalar self-coupling A^/^ and the top-quark Yukawa coupling ft- 
The gap equations yield a non- vanishing magnetic plasma mass for the 
gauge bosons, originating from the non-abelian self-interactions. We 
discuss in detail size and origin of higher order effects and conclude 
that the phase transition is weakly first-order up to Higgs masses of 
about 70 GeV, above which our calculation is no longer self-consistent. 
For larger Higgs masses even an approximation containing all g'^ con- 
tributions to Veff is not sufficient, at least a full calculation to order 
g^ is needed. These results turn out to be rather insensitive to the top- 
quark mass in the range mt = 100 — 180 GeV. Using Langer's theory 
of metastability we calculate the nucleation rate of critical droplets and 
discuss some aspects of the cosmological electroweak phase transition. 
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1 Introduction 



In the standard model the electroweak gauge symmetry is spontaneously bro- 
ken. However, at sufficiently high temperatures, above a critical temperature 
of about 100 GeV, the SU{2) x f/(l) symmetry of weak and electromagnetic 
interactions is restored Only many years after the ffist studies of 

symmetry restoration it was realized that in the standard model the rates of 
anomalous baryon- and lepton-number violating processes are unsuppressed 
at high temperatures This has important cosmological implications. In 
particular, it opens the possibility to understand, at least in principle, the 
generation of the baryon asymmetry of the universe within the standard 
model. Clearly, the study of "electroweak baryogenesis" requires a detailed 
knowledge of the transition from the symmetric to the broken phase in the 
standard model. Recently, this has led to a renewed interest in the elec- 
troweak phase transition 

Despite the fact that several important steps towards the understanding of 
the phase transition have already been made, we are still far from a complete 
description of this intriguing phenomenon. To a large extent this is due to 
the infrared problems of perturbative finite-temperature field theory. One 
manifestation of these infrared problems is the appearance of spurious terms 
in the effective potential Veff{(p), which are linear in the Higgs field if [0]. 
On the other hand, it has been argued (cf. [li H H) that the effective 
potential does not have a linear term. A related problem is the contribution 
of high order loop diagrams to low orders in the coupling constants. This 
can be understood in terms of dynamically generated plasma masses which 
damp infrared divergences. Much work has been devoted to the summation 
of daisy, superdaisy and other types of diagrams which yield higher order 



corrections to the effective potential 0-|TT|]. The free energy has also been 
evaluated in the 1 / A^-expansion and for background fields averaged over 
a finite volume [jl3[. Another important question concerns the decay of a 
metastable phase via nucleation, growth and coalescence of critical droplets. 
Furthermore, in the case of a weak ffist-order transition subcritical droplets 



and large thermal fluctuations are of importance [[I^ . 

In a recent paper |p, the phase transition in scalar electrodynamics was 
studied in detail. Using an improved perturbation theory, where plasma 
masses were incorporated from the beginning, the effective potential was 
calculated to order and A^/^, and it was explicitly demonstrated that the 
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spurious linear terms cancel. A complete set of one-loop gap equations was 
evaluated and used to determine a region in the plane of couplings (e^, A) 
where the symmetric phase is metastable. Following the theory of Langer 
the nucleation rate of critical droplets was calculated, and it was shown 
that a cosmological phase transition would have been first-order up to Higgs 
boson masses of the order of the vector boson mass. 

In this paper we extend the approach of ref. p to the standard model. We 
study an SU{2) gauge theory, i.e., the case g' = 0. This corresponds to the 
approximation where the M^-bosons and the Z-boson are degenerate in mass. 
We take into account the three generations of fermions, with left-handed 
doublets and right-handed singlets. The only relevant Yukawa coupling is the 
top-quark coupling ft. The effect of the other Yukawa-couplings is negligible. 

The outline of the paper is as follows. In sect. 2 we introduce our notation 
and study one- and two-loop contributions to the finite-temperature effective 
potential. We then discuss the breakdown of the ordinary perturbative ap- 
proach and stress the need for an improved perturbation theory. Sect. 3 is 
devoted to the structure of the vacuum polarisation tensor in spontaneously 
broken gauge theories at finite temperature. In sect. 4 we evaluate the ef- 
fective potential. We perform a ring summation, thus collecting all terms of 
order and A^^^, and we obtain a complete set of one- loop gap equations for 
the SU{2) gauge theory at finite temperature, which yield a non-vanishing 
magnetic plasma mass. Using the iterative solution of these gap equations 
we evaluate the effective potential up to cubic terms in the couplings. The 
equivalence of this method to the ring summation is shown, and some prop- 
erties of the potential are studied analytically. Sect. 5 contains a detailed 
analysis of higher order effects. We study the convergence of the improved 
perturbation theory, the influence of the obtained transverse plasma mass on 
the surface tension and the effect of superdaisy diagrams to order g"^ and A^. 
Following the theory of Langer [18| we compute the nucleation rate of critical 
droplets in sect. 6, where we also determine the nucleation temperature of 
a cosmological phase transition. A summary of our results and concluding 
remarks comprise sect. 7. 

Three appendices deal with matters peripheral to the main line of our 
discussion. In appendix A we prove that the ring summation with self- 
energy insertions at zero momentum yields all terms of order g^ and A^^^ in 
the effective potential. Appendix B gives the result of the "superdaisy" -type 
summation to order g^ and A^. Appendix C contains some details needed for 
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the computation of the nucleation rate. 



2 Conventional perturbative approach 

Let us consider the SU (2) gauge theory described by the Lagrangian C de- 
pending on a set of bosonic and fermionic fields, 

^ gauge ~l~ ^Higgs ~l~ ^fermion ~l~ ^ gauge fixing ~l~ ^ ghost- (l) 

Using standard conventions the bosonic part reads 

J^gauge + C^ggs = -\f;,F^^'' + {D,^)\D^^) - /.($t$) _ A(*t$)2^ 

/x<0, D^ = d^-ig''-^Wl a = 1,2, 3, (2) 
where the Higgs field $ is an SU{2) doublet. The fermionic part is given by 



jO-fermion = '^Li^D^IpL + ^ ^ j {iT2^*)tR + h.C. (3) 

Here r", a = 1, 2, 3, denote the Pauli matrices, ipL stands for all left-handed 
fermions and ft is the Yukawa coupling of the top-quark. All other Yukawa 
couplings are much smaller and can be neglected. The gauge-fixing and ghost 
Lagrangians read 

1 1 

C^gauge fixing = ^2^*^"^^" ' ^ ^^^^^ ~ l^VQ^Xa , (4) 

Cghost = CaM'^'c,, 5G'' = M-'5uj,; (5) 
here the SU (2) doublet $ of scalar fields has been written as 

where h is the Higgs field, (a = 1, 2, 3) are the three Goldstone bosons and 
9? is a constant background field. In the following we will use Landau gauge, 
i.e., we take the limit 77 — > 0. The classical minimum is at (f — v = ^J —n/X. 
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The tree-level vector boson mass m, the top-quark mass m^, the Higgs boson 
mass fh^ and the Goldstone boson mass rhy^ are given by 



m 



2 fV 2 f?,2 



4 ' '2 



< = (7) 



ml = Xi3^'-v% ml = Xi^'-v'). (8) 

Note that m| and are negative for small values of (p. The connection 
between couplings and zero-temperature masses reads g = 2mw/v, ft = 
^/2mt/v and A = mj^/2v'^, where v = 246 GeV. 

The finite temperature action is given by the integral 

Sf^= [dxC, (9) 

with 



/3 

f dx= f dr I d^x, P = ^. (10) 

Jl3 Jo J T 

Boson and fermion fields have to satisfy periodic and antiperiodic boundary 
conditions in r, respectively. 

The finite-temperature effective potential can be perturbatively evaluated 
in the loop expansion (cf. ref. [|l9l). The one- loop contribution is shown in 



fig. |l|. Its temperature dependent part can be calculated using standard 
techniques. The sum of the leading term in the high temperature expansion 
and the tree-level potential yields: 



VM + VI'\^,T) = i '_1±^^T^ _ XvA + '-X^\ (11 




As is well known, this potential predicts a second-order phase transition. 
Note, that the finite-temperature contribution due to scalar loops is propor- 
tional to A and therefore of the same order as the tree-level term A(y9'*/4. 
Going one order further in the high-temperature expansion one obtains 

T/r ^^T/(2)r 1 /3^^+8A + 4ii 2 , 2\ 2^1, 4 
K)(^) + ^1 (^,r) = 2 I 16 -At; +4^*^ 

- (9m3 + ml + 3mj) ^ . (12) 
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The potentials (|TT]) and ( [I^ ) differ in several important respects: 
(i) We are faced with the familiar problem that fh^p and become imaginary 
for small values of (p, i.e., the naive one- loop effective potential is complex. 
This phenomenon is not an artefact of the high-temperature expansion, since 
the exact one-loop contribution of the Higgs field is proportional to 

J+(mJ) = dxxHog (l - e-V-^+(3^^-^)VT^) ^ (13) 

which is also complex for small values of y?. A similar expression is obtained 
for the Goldstone fields, and no cancellation appears between the two contri- 
butions. At zero temperature, the imaginary part of the potential has been 
related to the lifetime of a particular quantum state ||20l . 



(ii) The new cubic terms are of order g'^ and A^''^, and thus of higher order 
than the finite-temperature corrections in Vi^\ Clearly, these new contribu- 
tions must be combined with terms of the same order which appear in higher 
order loop diagrams. 

(iii) The cubic terms originate from integrals which are infrared divergent 
for vanishing masses. No cubic term results from the top-quark loop, since 
the fermionic modes always have non-zero Matsubara frequencies {un = 
(2n + l)7rT), and hence do not suffer from infrared problems. 

(iv) In order to see the qualitative features of the effective potential (0) 
we first ignore the terms proportional to A^/^, which can be justified for 
A <^ g^. For temperatures above a certain "barrier" temperature T^, the sec- 
ond derivative of this modified potential is positive, thus the symmetric phase 
is a local minimum of the effective potential. Here the barrier temperature 
Tb is 

In eq. (|12D the terms proportional to (f^ and ip'^ have opposite signs. As 
a result of a compensation between these two terms, for temperatures close 
enough to Tf,, another minimum exists which, at a critical temperature T^, 
is degenerate with the minimum at (p = 0. Hence, in this approximation the 
phase transition is first-order. 

According to (ii) we have to consider higher loop terms involving boson 
fields, first the two-loop contributions shown in fig. ^ These graphs yield 
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linear terms in the temperature dependent part of the effective potential: 
yiscalar)^^^^^ = - ^ (3^ + 8 A + 4/,^) (m, + SmjT^ 

+0{g\X') , (15) 

vt''"'''\v,T) = -^g'^T' + 0{g\X') . (16) 

oZn 

The expressions are cubic in the couplings, and therefore they must be com- 
bined with the contributions of the same order in Vi {(p,T). However, the 
combined result does not solve any of our previous problems related to the 
scalar sector, but rather introduces a new problem, the appearance of a very 
disturbing linear term in the effective potential. The fact that the two-loop 
diagrams yield terms of order g^ and X^^"^, and not only g'^ and is another 
manifestation of the infrared problems in finite-temperature field theories. 
The two-loop contributions can be considered as one-loop graphs where the 
other loop plays the role of a self-energy insertion. The infrared divergencies 
are cut off by these dynamically generated masses (proportional to g or a/A), 
which leads to factors T/rrii and thereby reduces the order in (7 or A of the 
corresponding diagram. 

In appendix A we show that to order g^ and A^^^ only ring diagrams con- 
tribute (cf. fig. 0). The complete summation of these terms with the proper 
combinatoric factors is required to obtain the correct effective potential to 
order g^ and A'^^^. Therefore, we need the self-energy terms at zero exter- 
nal momentum. The evaluation of these contributions is straightforward for 
scalar fields, but for vector fields the polarisation tensor needs special atten- 
tion in spontaneously broken gauge theories. The next section is addressed 
to this question. 



3 The structure of the gauge boson propaga- 
tor 

Most of the results of this section have already been discussed in ref. [|]. 
Nevertheless, they are included here to make our paper self-contained. 

In order to determine the plasma masses of vector and scalar fields, we 
first have to discuss the structure of the vector propagator at finite temper- 
ature. The gauge boson self-energy n^,y(fc) depends on the 4-momentum k'^ 
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and the 4-vector = (1,0) which specifies the rest frame of the system (cf. 
ref. [0). Hence, in general II^j, is a hnear combination of four tensors. A 
convenient choice is 

Pt,. = 9;U,-^]9\, (17) 



PL^,u = '^^-g^u-PT^,u = pK^l^ (18) 
k,, k,. 



Po.^ = (19) 



7^ (^ul + k,ul) , (20) 



where = — k^^ is transverse, uj^k^ = 0. These tensors satisfy the 
relations 

P^ = -Pt, PI = -Pl, Ph = -PG, = ]^{Pl + Pg),{2^) 
PtPl = PtPg = PlPg = SPt = PlSPl = 0, (22) 
Pt/ = 2Pl/ = 2Pg/ = -2, 5'/ = = . (23) 

In Landau gauge, where the bare propagator is transverse, 

Dik) = ^^,{PL + PT), (24) 
the gauge-boson self-energy tensor 

U{k) = UL{k)PL + UT{k)PT + Ils{k)S + IlG{k)PG (25) 
yields the full propagator 

oo 

D{k) = D{k)mk)D{k)r 

n=0 

-Pl + n.. PT- (26) 



k^ -m^ -UL{k) " k^ -m^ -Uriky 



Due to the relations (piD and (|22D the full propagator does not depend on 
Usik) and Udk). 
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However, knowledge of Hg is important since it enters in the relations 
which yield the longitudinal and transverse plasma masses: 

6ml = Ul{0) = Tr[n(0)Pi] = -noo(O), (27) 
6ml = Ut{0) = ^Tt[U{0)Pt] 

= -l[n%(o) + ni(o) + nG(o)] , (28) 

where 

Ucik) = Tr[n(A;)PG]. (29) 
Here ll^^{0) is defined by setting first ko = and then performing the limit 
A;^ — i> 0. In gauge theories with unbroken symmetry one has H^ = in Lan- 
dau gauge. Note, that this is not the case if the symmetry is spontaneously 
broken. As eq. (^) shows this fact, which seems to have gone unnoticed in 
the literature, is important in order to extract the correct transverse mass 
from Hf^i^. As we will see, in non-abelian gauge theories this transverse mass 
turns out to contain a field independent magnetic plasma mass which plays 
an important role in the estimation of the higher order effects. Hence, the 
unambiguous determination of this term is needed. 

4 The effective potential in next-to-leading 
order 

A. Ring summation 

The one-loop gauge boson self-energy corrections for the SU (2) gauge theory 
are shown in fig. (^. The self-energy contributions have very complicated 
dependence on the momentum (fco, |^|)- As we show in appendix A, to order 
and X^^^ one has to set fco = for the external lines of the self-energy 
diagrams in the ring summation. On the other hand, by use of a Taylor 
expansion in \k\ it is easy to show that only the \k\ limit of the self- 
energy graphs contribute to the effective potential in this order. In this limit 
we find for the longitudinal and transverse plasma masses to order and 

6ml = ^9^T^-^(^^ + m^ + Sm^ + lQm]T, (30) 
6 loTT \m + m^ ^ A. J 
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Here we have neglected terms of higher order in the high-temperature expan- 
sion, since we will not need the corresponding terms in the effective potential. 
It is worth mentioning that no field independent transverse plasma mass ap- 
pears in eq. (|3TD, in accord with the well-known fact that in one- loop order 
it vanishes in all gauges. Having also determined the corresponding scalar 
self-energy contributions, one can perform the ring summation. The corre- 
sponding diagrams are shown in fig. ^ where the "blobs" stand for one- loop 
self-energy contributions (cf. fig. ^and similar scalar terms). Special care is 
needed with respect to the infrared limit and the combinatoric factors, partic- 
ularly in the case of one self-energy insertion, since these graphs are two-loop 
diagrams with different symmetries depending on the type of propagators. 
The ring summation yields the potential 

-(3mi + 6m| + mj + SmJ)^ + 0{g\ \\ ff), (32) 

where is given by eq. (0) and the masses are the sum of tree-level terms 
and one-loop self-energy corrections to order g"^ and A 

2 11 2T.2 , 

= -gT 

2 

< = (^^' + ^ + i/')(^'-^')+V. (33) 

The potential ( |5^ describes a first-order phase transition, since there ex- 
ist two degenerate minima at a critical temperature T^, similar to the case 
considered in sect. 2. However, due to the non- vanishing longitudinal plasma 
mass the strength of the transition is less than the high-temperature expan- 
sion of the one-loop result suggests. Clearly, a non-zero field independent 
magnetic mass would further weaken the transition. 
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The potential ( ]32D contains no term linear in Hence, the linear terms 
which appear in the two-loop expression are cancelled by contributions from 
the ring summation. As we show in appendix A the result contains all 
terms of order and A^^^. Therefore the appearance of linear terms of order 

or 

A3/2 

can be ruled out. To higher order in g and A the absence of a 
linear term has not yet been proven. Note, that for temperatures above the 
barrier temperature Th all masses, and therefore also the potential, are real. 

Formal arguments concerning the existence or nonexistence of linear terms 
in the framework of the perturbative expansion in the coupling constant are 
questionable, since such an expansion does not necessarily reflect the be- 
haviour of the potential at the origin. As an illustration consider a hypothet- 
ical potential containing a term proportional to 

^g2n+2rp2 + g2n^2_ (34) 

This function is symmetric in (f and contains no linear term. Nevertheless, 
a formal expansion in powers of g gives rise to a linear term, 

n+2rp2 

Jg2n+2rp2 + g2n^2 ^ gU^ ^ y__^ ^ (35) 

On the other hand, a potential with a term 

^^g2n+2rp2 + g2n^2 (35) 

does contain a linear term which, however, is not visible in the expansion in 
powers of g: 

^^g2n+2rj^2 + g2n^2 ^ gn^2 ^ ^"+2^2/2 + C(^"+^). (37) 

The behaviour of the potential obtained by an expansion in the coupling 
constant does not coincide with that of the exact potential. Note, however, 
that in both cases the ratio of the next-to-leading term and the leading term 
becomes very large at ~ which indicates that the perturbative expansion 
in g is not valid in this region. Because of this fact we will carefully study 
the effects of the higher order terms in the couplings at ^9 ^ and close to 
the second non-trivial minimum of the effective potential (^2]). 

The ring summation yields a sensible result to order g^ and A^^^. This is 
rather surprising since individual terms in this sum are not well defined for 
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small values of due to the imaginary masses which appear in the scalar field 
propagators. This inconsistency would show up in the next higher order, i.e., 
g'^ and A^, where all diagrams have to be added to the ring diagrams which 
contribute to this order. 

B. Gap equations 

In this subsection we shall evaluate the potential ( ^2|) in a more elegant 
way which will also enable us to estimate higher order corrections to our 
result. For small values of the temperature dependent plasma mass cor- 
rections 5m\ rp can become larger than the tree-level mass m^. This suggests 



an improved perturbation theory g, |2T[ , where loop diagrams are evaluated 
with boson propagators containing the exact masses m| = + 5m\ j,^ 
m'^^^ = fn'^^^ + Sm'^ .^. The radiative corrections are treated as counter terms 



SSf3 = -^J^dp[A^{p){6mlPL^. + 6mlPT^.)A^{p) 

-6ml^\p)-6mlx\p)], (38) 



and are determined self-consistently by solving gap equations at the corre- 
sponding loop order. As it was pointed out above, there is no need to use 
this procedure for fermions, since fermions do not suffer from infrared prob- 
lems. The ghost propagator remains massless in Landau gauge. Therefore, 
in our approach the bare fermion and ghost propagators play the role of full 
propagators. The one-loop contributions to the effective potential with these 
dressed propagators are shown in fig. |^. Note, that here and in the follow- 
ing the blobs stand for exact masses and not only for one-loop self-energy 
insertions. The mass counter terms for top-quark and ghost are zero. 

The relevant one-loop graphs for the gap equations are shown in fig. ^ 
A rather lengthy, but straightforward calculation yields the following set of 
equations to order and X^^"^: 

ml = ^g^T^ + ^^-T^ + + 3m^ + IQmA T, (39) 

2 g^T . g^ ( 

mj, = — — mx + m — — 



Qtt \mT + m^p 




2 _ / I I ^ £2] (rr2 rr2\ , Q-2 



16 2 4" 



12 



3^7 



2 r 



16n 
3A 

An 

167r 



2/1 2 , 

+ 2mT + m ( 1 



T 



2/ 3 1 
+ + m { \ 



T, 



(41) 



+ ^ + If^) - T^) + 
{rriL + 2mT) T 



m 



T, 



where 



m 



9^/2, 



m 



(42) 



(43) 



' 3g^ + 8X + Afr 

As aheady mentioned, the ghost mass remains zero, and the mass corrections 
for fermions are not important since their Matsubara frequencies are always 
at least (9(1) -T. Thus, there is no need to consider gap equations for fermions 
and ghosts. 

It is instructive, first to consider the gap equations for the pure scalar 



theory {g = 0, /t = 0) at = 0. From eqs. (|l|) and (|42|) one obtains for 
temperatures close to the barrier temperature 



m<^ ~ ~ (T - Tb 



(44) 



i.e., and approach zero with critical index one. This well-known result 
was first obtained by Dolan and Jackiw in the large- iV limit 0. We obtain 
the same result, because at ip = only graphs e) and /) of fig. | contribute, 
which are the leading terms in the 1/A^-expansion. 

Of particular interest is eq. (^Oj) for the magnetic mass. At ip = 0, one 



has = m^, and therefore 



my 



9^^ 
3n 



■rriT- 



(45) 



This equation has two solutions. 



rriT = 0, 



(46) 
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and 

mr = ^, (47) 

thus an unwanted ambiguity seems to appear. However, only the second 
solution is physical, since only = g'^T/{37r) can be continuously connected 
to a positive solution of eq. ai ip > 0. The (p limit of the negative 
solution of eq. (^) is ttit = 0. This negative solution is unphysical, since rriT 
is defined as absolute value of a square root appearing in the high temperature 
expansion. The result (|4^) has been independently derived in ref. p. 



The obtained value rriT = g'^T/{37r) for the transverse plasma mass goes 
beyond the accuracy of our calculation which is only valid to order g. How- 
ever, the appearance of a magnetic mass of order g^ is expected to be a 
nonperturbative feature of the theory ||22|. In order to study the depen- 



dence of our results on the unknown value of the magnetic mass we will take 
rriT = 7(7^T/(37r) as the ip —>■ limit. 

It is well known that in the abelian Higgs model the transverse plasma 
mass is zero at finite temperature, which is in agreement with our previous 
results 0. In nonabelian gauge theories the situation is different. In this 
case the gauge boson self-couplings (cf. figs, ^q, ^) yield a non-vanishing 
transverse plasma mass. This magnetic mass can be self-consistently deter- 
mined by summing an infinite set of diagrams. The non-zero solution (^) for 
the transverse plasma mass can be viewed as a result of an infinite iteration 
process with an arbitrarily small initial value. In a diagrammatic picture 
the iterative solution corresponds to an infinite sum over "superdaisy"-type 
self-energy insertions. The appearance of a magnetic mass term has been 
studied by other methods in finite-temperature QCD ||2^. Lattice SU{2) 
models and infinite summations also give values for 7 which are 0{1). 

Let us now solve the gap equations, which represent a set of nonlinear 
equations with four variables: 

M'^ = f{M), M = {mL,mT,m^,m^) , = (m^, m|, mj, mj) . (48) 

Except for linear problems, finding of roots invariably proceeds by iteration. 
Starting from some approximate solution, a useful algorithm will improve its 
accuracy. We are interested in a Taylor-expansion in g and of the result. 
Setting g = = gives zero for Ai, thus the Taylor expansion for Ad starts 
with g or -\/A, and for Ai^ with g'^ or A. Suppose that we know the results for 
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A4 to order n. Inserting this trial solution in the right-hand-side of the gap 
equations we obtain an expression to order n + 2 for M"^, and to order n + 1 
for Ai. The coefficients of the previous orders do not change. Consequently, 
the above iterative procedure reproduces the Taylor-expansion of the exact 
solution order by order. 

This iterative solution has a diagrammatical picture. The lowest order 
result corresponds to the one-loop self-energy insertion at zero external mo- 
mentum. A new iteration means one order higher in the coupling, e.g., the 
sum of the one-loop self-energy insertions of the one-loop self-energy inser- 
tion, the so-called "daisy" graphs, the next order gives a new self-energy 
insertion on the bare lines, and so on and so forth. After an infinite number 
of iterations one obtains the sum of all "superdaisy" diagrams [0]. Note, 
however, that this summation does not give the correct answer for the ef- 
fective potential to order g^, and higher. One cannot justify the zero 
external momentum limit of the gap equations for higher order contributions 
(cf. appendix A). 

Our goal is to calculate the effective potential to order and A^/^. Due 
to the global SU{2) symmetry of the theory the potential is only a function 



of v2$t$ = y't^z _|_ ^2 _|_ ^2 _|_ fjgnce, at X = 0, the masses m^{(p, T) and 
m^{ip,T) are given by 



m>,T) = -^^. (50) 

Thus, in order to obtain the potential to order and A'^''^ one has to evaluate 
or to the same order. This means two successive self-energy insertions, 
which is exactly the ring summation. The result reads: 

Vgap{ip,T) = / dip' Lp'ml{(fi',T) 



2 V 16 2 4''*; ^ wv- ^ 

+0{g\\\ft), (51) 
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where the lowest order masses are 



(1)^ 11 2rr2 , 2 

my = — q I + m , 

(1)2 2 1 4rr2 I 2 

"^T = 7 ^fi- T +m , 

= i^Q9' + ^ + \miT'-T^) + 3m\ 

mf = {^/ + ^ + \j?){T'-T,')+m\ (52) 

As it should be, this result is identical with the ring potential (|32D for 7 = 0. 
In order to illustrate the dependence of our results on the unknown value 
of the magnetic mass we have modified the iterative solution of the gap 
equations to leading order in the couplings in such a way that rrij^^ {(p = 0) = 

The evaluation of the effective potential to order g'^, and requires 
the incorporation of two- loop contributions in eqs. (^) - (^2]), and in order to 
obtain the exact effective potential one has to solve the full Dyson-Schwinger 
equations, a non-trivial task! However, the obtained gap equations already 
contain some part of the higher order terms, namely the "superdaisy" dia- 
grams, and thus contain valuable information about higher order corrections. 
Before we investigate these questions in the next section, we study some prop- 
erties of the obtained potential (|5T|) . 

C. Plasma masses and the order of the phase transition 

A nice feature of the one- loop effective potential considered in sect. 2, 
with or without cubic term, was the possibility to study the phase struc- 
ture analytically. Due to the plasma masses our potential (|32D or (^TJ) now 
contains terms of the form (C^ + yj^)^/^, and is therefore more complicated. 
Nevertheless, it is still possible to study its main properties analytically. 

In order to understand the effect of plasma masses consider a hypothetical 
potential with one plasma mass term proportional to T, 

V^T, V) = ^iT^- T^)V' - y(^'T^ + V"f'" + I^'' (53) 

where a, b, c and A are positive real numbers, independent of temperature. 
Vhyp has a local minimum at the origin as long as its second derivative there 
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is positive. This condition yields a corrected barrier temperature Tb, 

The phase transition described by this potential is first-order if at the 
barrier temperature Tf, a second minimum exists at non-zero (p, otherwise it 
is second-order. The condition that at the first derivative of the potential 
with respect to (p vanishes at some value ip > implies for the couplings 

A<|. (55) 

In the standard model only the parameter a depends on the top-quark 
Yukawa coupling, hence this condition does not depend on the top mass. 
Clearly, it implies that at a critical value of the Higgs mass the transition 
changes from first-order to second-order. It is easy to show that for several 
plasma mass terms of the form 

^(c?T^ + ^y/' + ^-ficlT' + p?f^ + ... + ^(c^T^ + ^^)^/^ (56) 



the condition (|55| ) for a first-order phase transition becomes 



A<^ + ^ + ... + ^. (57) 

2ci 2c2 2c„ ^ ' 

If one neglects in the electroweak potential (|3^ ) terms of order A^/^, thus 
assuming ^ A, the cubic term with longitudinal plasma mass is sufficient 
to give a first-order transition if vnn < 13 GeV. However, because of the 
second cubic term with zero plasma mass the phase transition is first-order 
for all Higgs masses. This situation changes if a transverse plasma mass is 
included, as in eq. (pi]). The corresponding term in the effective potential 
yields a second-order transition if 

85 GeV , , 

TTlH > ^ . (58) 

Inclusion of the cubic term with longitudinal plasma mass slightly relaxes 
this bound according to condition (|57D. As already mentioned, 7 = 0(1). 
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The scalar part of the electroweak potential needs a different treatment. 
Set g = 0. Now the potential is of the form 



+ {a{T^ - T^) + 3A(^2)3/^)t, (59) 
and the temperature, where the potential barrier vanishes, is 

^' - T^WxTa 

The only solutions of the equation dVhyp{Th, {p)/dip = are 

<^' = 0, <^2 = 6AT,2(-4±z), (61) 

i.e., the scalar part of our potential always yields a second-order phase tran- 
sition, independent of the value of A. 

In the general case of non-zero plasma masses for vector and scalar bosons 
(cf. eq. (^Tj)), the order of the phase transition changes at a critical Higgs 
mass which can be evaluated as function of 7. The result is shown in fig. 
(0). Note, that for the present lower experimental bound on the Higgs mass 
the phase transition is always second-order for 7 > 2.5. 

As we shall discuss in detail in the next section, also the improved per- 
turbative approach becomes inapplicable close to the barrier temperature. 
Therefore, we cannot calculate the effective potential at this temperature. 
Nevertheless, many features of our numerical analysis will be very similar to 
the result of our analytical study in this section, which one may expect based 
on continuity arguments. 



5 Higher order effects 

In this section we will estimate the size of different higher order corrections 
to the effective potential (pi]). First we will use an extended version of the 
method discussed in ref. @ to check the convergence of the perturbation 
series. We will then study effects of a non-vanishing magnetic mass, in par- 
ticular its influence on the surface tension, and finally we will perform an 
analysis including some higher order terms in the effective potential. 
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A. Convergence of the perturbation series 



At the origin, ip = 0, the curvature of the potential m^(0,T) reads (cf. 
eq.Q) 




3 2 A 1 

—9 H h - 

16 2 4 



:TF9' + 7; + ^miT'-T^) , (62) 



27r V 16 2 4^^ 

and the condition of vanishing curvature, m^(0,T^') = 0, yields the corrected 
barrier temperature 



rr = + 8A + Af! - + 0{g\ A^ //) j . (63) 

Note, that there is no correction of order A^^^. 

At ip > 0, the gap equations ( pQ]) - (|i2|) contain terms proportional to 



T/rrii which reflect the expected infrared problems of flnite-temperature per- 
turbation theory. The perturbative expansion is reliable if these terms are 
smaller than leading-order terms. Inspection of eqs. ([39| ) - (^2]) shows that 



this is guaranteed if the following inequalities are satisfled: 

(64) 




- + -)<! ■ (65) 

47r \m^p 

The first inequality stems from the equation for m\ and the second from 
the equation for m^. We have included a factor ^ which ensures that leading 
terms of order g"^ and A are ^ times larger than next-to-leading contributions. 
Clearly, for fixed values of g,X,ft,T and ip, the larger the value of ^, the 
better the behaviour of the perturbative expansion. If the conditions (|6^) 
and (|65| ) are satisfied one expects the uncertainty of the field dependent 
part of the effective potential ( pID to be of order 1/^^. We will now study 
the implications of the above conditions in two ways. First we will use the 
iterative solution of the gap equations and insert the squared plasma masses 
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to order and A (cf. eq. (|^)) into the conditions (0) and (|65D. We will 
then repeat the analysis in a more conservative way by using the squared 
plasma masses in the ^-conditions to order and A^/^, which are smaller 
than the ones to leading order (cf. eq. (|49|)). 

Close to the origin, at ~ 0, the above conditions imply that one cannot 
even reach the barrier temperature T^. For 7 = eqs. (|52D, (|64D and (|65|) 
yield the lower bounds on the temperature T, 

T>Ty , T>Ts , (66) 

where 

~ = ^M. (671 

and 

rp 1^2 C2\2 

Ts' 3vr2(A^2 + | + l/,2)- ^''^ 

Here, the subscript "l^" ( "S"' ) indicates that the infrared divergence for the 
vector (scalar) field plasma mass sets the temperature. Hence, our expression 
for the effective potential is only reliable for temperatures above T*, which 
denotes the largest temperature among Tl,Tv and T5. Non-zero values of 
7 give a non-zero magnetic mass rriT at ~ 0. The corresponding infrared 
behaviour leads to a smaller temperature Ty. 

In order to establish the existence of a first-order phase transition for a 
given value of the scalar self-coupling A we have to show that for some values 
(p > 0, where the ^-conditions (|64D and (|65|) are satisfied, 

V{ip,T*) <V{0,T*) . (69) 

Note, that T* = T*{g,\), and thus the ^-conditions are satisfied at ip = 0. 
The larger the Higgs mass, i.e., the scalar self-coupling A, the more difficult it 
is to fulfill the condition for ^5. Let us now consider the potential in a region 
> 0, where the .^-conditions are satisfied. At a certain maximal value of A, 
and a corresponding critical temperature T*, the minimum of the potential 
in this region, at > 0, is degenerate with the minimum at = 0, 

y(^.,T;) = \/(o,T;). (70) 

Note, that in general T*{g, A) < Tc{g, A), the true critical temperature, and 
Vfig, A) > ipcig, A), the local minimum at Tc. Hence, ^p^, is not necessarily an 
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extremum of the potential with respect to the full range in ip. As discussed 
above, the condition T* = T* defines a maximal value of A for a given value 
of ^. We have plotted the corresponding functions rriH^^s) and mfji^v) in 
fig. p. The condition obtained from the gauge boson gap equation is always 
satisfied with > 7.6, even for 7 = 0. The result is essentially independent 
of the top mass in the range rrit = 110 — 180 GeV. Both, and are larger 
than 2 up to Higgs masses of approximately 200 GeV. 

Let us now repeat this analysis using the higher order plasma masses (|^) 
in the ^-conditions. The allowed minimal temperature T* is higher than in 
the previous case. Since the Goldstone mass = 2dV/d(p^ vanishes at the 



second minimum Lp(T*), the condition (|65|) cannot be satisfied here and the 
perturbative approach breaks down. Hence, one always has > ip(T*). We 
can now determine the temperatures T*{g, A). In order to illustrate how far 
T* lies below the true critical temperature T^, we have plotted an effective 
potential at T* in fig. |. The condition T*{g, A) = T*{g, A) again defines a 
line in the ^-A-plane, denoted by in fig. |^. The difference between the 
two boundaries and in fig. | is considerable. For Higgs masses above 
80 GeV the value of is smaller than 2. 

B. The effect of a magnetic mass 

Our estimate of the largest Higgs mass up to which perturbation theory 
is reliable, which we have carried out in the previous subchapter, depends 
on the parameter ^. The appropriate value of ^ can only be determined by 
calculating higher order corrections. In sect. 4C we saw that the magnetic 
mass has an important effect on the phase transition. Hence, the effective 
potential has to be calculated at least up to terms proportional to g^, the 
order at which the magnetic mass contributes. Since we have not calculated 
all graphs contributing to this order, we have introduced a constant 7 in eqs. 
( p2D which parameterizes the uncertainty in the size of the magnetic mass. 
The solution of our gap equations gives 7 = 1, other approaches (cf. ref. 
[^) yield values of the same order of magnitude. 

In order to estimate the effect of the magnetic mass on the phase transi- 
tion we calculate the surface tension, a physical quantity, from the potential 
( PI ) as function of A and 7: 

^{Kl) = j^^ d^pV{^,%), = if{T,). (71) 
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As discussed in the previous section, the temperature Tc is always larger than 
the temperature T*, if the chosen value of A is allowed by the value chosen 
for ^. Therefore one can reliably calculate the potential close to the origin, 
(fi ^ 0, and for large values of (p. Of course, in the intermediate regime, which 
is needed in order to get the surface tension from eq. ([7l|) , the potential has 
large uncertainties. Nevertheless, as the dominant contribution comes from 
the vector loops which have no infrared problems in the intermediate range of 
(fi, it is conceivable that the surface tension, which is an integrated quantity, 
can be obtained from eqs. (^) and (|7TD to good approximation. 



We have plotted a as function of the Higgs mass for different values of 7 in 



fig. |10|. Note, that the disappearance of the surface tension means transition 
to a second-order phase transition. The Higgs masses for which a vanishes 
for different values of 7 are in agreement with fig. ^ For a given value of tuh 
the effect of the magnetic mass on the phase transition can be characterised 
by the ratio 

C(7) = 1 7 ^ — 7 ^ • (72) 

Clearly, if (^(7) <^ 1 unknown higher order corrections are so large that a 
first-order phase transition cannot be established unequivocally. According 
to fig. |10] the quantity ( is of 0{1) for Higgs masses up to about 70 GeV. 
This value of mn is rather close to the present lower experimental bound. 

C. Scalar loop and superdaisy contributions 

In this subchapter we analyse the role of some higher order terms in the 
perturbative expansion. As in the previous subchapter we again study the 
surface tension ([TI]). 

It is instructive to first set g = and to look at the perturbative expansion 
in A. The phase transition of the pure scalar theory is well known to be 
second-order. Our results agree with this fact, since the effective potentials in 
the leading order, C(A), and the next-to-leading order, 0(A^/^), give second- 
order phase transitions. The appearance of a first-order phase transition 
would only be a sign that our perturbation expansion is not reliable, since 
it would result from a cancellation between terms of different orders. In 
other words, contributions from different orders would be of the same order 
of magnitude. 
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Let us now study the case 7^ 0. The scalar contributions together with 
the term 0{g^) give a second-order phase transition, while inclusion of terms 
of order yields a first-order transition. The above mentioned problem is not 
relevant here, since the cancellation takes place between terms proportional 
to different expansion parameters, namely and A. In this perturbative 
approach, where g^ ~ A, the size of the next-order terms (in our case A^^^ 
and g^) must be small, and in particular their effect on the surface tension 
has to be a reasonable correction. 

We have determined the contribution 0(A^/^) to the effective potential. 
In fig. ^ we have plotted the surface tension at the critical temperature 
as function of the Higgs mass for the potential containing terms 0{\^^'^, g^) 
(solid line) and for the potential with terms 0{X,g^) (long dashed line). For 
small Higgs masses the difference is small, but at approximately 100 GeV 
the difference becomes 0(50%) indicating a breakdown of the perturbative 
approach. 

Superdaisy diagrams can be summed by finding exact solutions to the gap 
equations. Iterating the gap equations we have determined the contribution 
of the superdaisy diagrams to order g'^ and A^. The rather lengthy expres- 
sions for the potential are given in appendix B. In fig. 11 we have plotted 
the surface tension at the obtained critical temperature for this effective po- 
tential (short dashed line). The superdaisy terms give a correction 0(20%) 
correction, and thus do not change the qualitative picture of a first-order 
phase transition. 



6 Decay of metastable states 

In condensed matter physics the decay of metastable states is described by 



hanger's theory [|T8|. The starting point is a coarse-grained free energy which 
depends on the order parameter, temperature and the coarse graining scale. 
As function of the order parameter the free energy has a local metastable 
minimum which is separated from the global minimum by a barrier whose 
height determines the lifetime of the metastable state. From a stationary so- 
lution of the Fokker-Planck equation for the probability distribution of large 
fluctuations of the order parameter ( "subcritical droplets" ) one then obtains 



a formula for the decay rate [18| which depends on the free energy of the 



metastable state, the free energy of a saddle point field configuration which 
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interpolates between local and global minimum and a "dynamical factor" 
which cannot be obtained from equilibrium thermodynamics. 

Langer's formalism can be directly applied to the decay of metastable 
states in quantum field theories. The case of scalar electrodynamics has 
been studied in a recent paper p. The decay rate is given by 

^_ K ImZf^m ^^^^ 



2tx Zfi{^ = 0] ' 
where 

Z,l<S>] = e-/5^l*.r] (74) 



The integration measure for the vector field includes the gauge fixing and 
ghost terms, and stands for all fermion fields. In general, the free energy 
also depends on and ip. However, we will only consider stationary points 
with Afj^ = ip = 0. For simplicity we have therefore omitted the dependence 
of F[$,T] on Afj_ and as well as the corresponding functional derivative 



terms in the exponent of eq. ([7^) . $ is a field configuration which inter- 
polates between the symmetric and the broken phase. Since $ = and $ 
are approximate stationary points of the free energy F[^,T], we neglect the 
second term of the integrand in eq. (ffil). 

The functional integral over vector and fermion fields yields an effective 
action which depends on the scalar field $ + $, 

exp ( - / dx(dj<^ + <l>)^9^(<l> + <!>)- Vo(z) 

Jf3 

- V{z, T) + Z{z, T)9^($ + ^)^d^{^ + <!>) + .. .)), (75) 



/3 



where z = y2{^ + $)t(<|> + $). Note, that this expression is invariant under 
the global symmetry 0(4) rather than SU (2) xU{l) which is the symmetry of 
the Lagrangian given in eq. (|l]). Indeed, integrating out the top-quark yields 
additional wave function correction terms Zi{^ + ^) which are invariant only 
under the smaller symmetry group. For simplicity we will neglect all wave 
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function correction terms in the following. A more detailed discussion will 
be given elsewhere 1 2^ . 

At one-loop order the potential V{z, T) is well known. Including vector 
boson and top quark loops one obtains for the full potential in the high 
temperature expansion (cf. eq.( 



32n 



(76) 



The potential has two local minima, at 93 = and at ^p{T) > 0. At the 
critical temperature Tc both minima are degenerate and one has 

= ^{Tc) = ^^^c- (77) 



The integral over the scalar field fluctuations $ can now be carried out 
in the saddle point approximation. In the thin wall approximation ||2^ the 
stationary point of the approximate free energy 



F[$,T]= d^x{\V^\^ + V{z,T)) 



(78) 



which appears in the integrand of eq. (|7^ , has a saddle point $ which 
interpolates between $ = and <I>(T) > 0. For temperatures just below the 
critical temperature the height of the barrier between the two minima is 
large compared to the potential difference between the minima. In this case 
the saddle point can be computed in the thin wall approximation and one 
obtains 



<l>(r) 



with 



1 — tanh ( - 



d 



(79) 



a 



dip 



R{T) 



2a 



e{T) = V{0,T)-V{^{T),T). 



(80) 
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The free energy Ftw[^, T] in this approximation is then the sum of a volume 
term and a surface term: 

- - Att 

Ftw[^,T] = 4rf(T)(7 - —R%T)e{T). (81) 

It is sufficient to evaluate the surface tension a at the critical temperature 
T 

For the partition function of the saddle point one now obtains 

Im = ^^vn^(A,)-V2e-/^^[*'2^], (82) 

A_ i 



where denotes the product of all positive eigenvalues Aj of fluctuations 
around A_ is the single negative eigenvalue, and V is the volume of zero 
modes associated with the symmetries of the system under consideration. 

The scalar fluctuations $ consist of the radial modes (f and the Goldstone 
modes xf. 

F[$ + I>,r] = F[^,T] 

+ i ^ dx (fi{x){-A + U^{r))<f{x) 

+^ dx Xi(x)(-A + U^{r))Uy)- (83) 

The corresponding potentials for the scalar flelds (p and % are 

Ux(r) = -Anv,T) (86) 



ip dip 



The spectrum of eigenvalues contains six zero modes, three for translational 
invariancc and three for the global SU (2) x f/(l) symmetry of T] which 
is spontaneously broken to the electromagnetic U{1) subgroup. The cor- 
responding volume factor for the translational modes is well known. The 
volume factor for the global symmetry can be calculated in the same way as 
for sphaleron tunneling processes, the details are given in appendix C. One 
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obtains: 



V, = tiizld'^A"'^ (86) 



2 V27r 

where V is the total volume of the physical three dimensional space 



The discrete <^ spectrum is well known ||18|, since the bound states 



are localized at r i?. There is one negative eigenvalue, 

A- « - Jj, (87) 

which guarantees that ^/3[$] is purely imaginary. Furthermore, there are 
"Goldstone modes" which correspond to deformations of the droplet surface 



[p!q , [25| , [26[| . The corresponding contribution to the determinant of eigenval- 
ues is where ^i = m^(0,T) 0. Combining eqs. (|73|), (||) 
and (0) we finally arrive at the transition rate 

V - 2OT$ ^ '"-'^'^ (^"^"'^ («^'"'" 

Here the contributions of zero modes and Goldstone modes to the deter- 
minant of scalar fluctuations around the saddle point have been taken into 



account. The "dynamical factor" k has recently been evaluated 

IGrjcr 



(89) 



where Auj is the difference of the enthalpy uj = —TdV/ dT between symmetric 
and broken phase, and the viscosity rj = 65.4 in the standard model . 
The pre-factor of the exponential in eq. (|88D and the naive estimate 
turn out to give in our numerical calculations essentially the same results. 
However, a priori this is not clear, and the evaluation of the pre-factor is 
needed to verify the validity of the semiclassical approximation. 

Eq. (|88|) gives the decay rate of the metastable symmetric phase in the 
framework of hanger's theory of metastability. The free energy F[$,T], 
obtained by integrating out vector boson and fermion fields, plays the role of 
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V 


V2/3 


Vx 


V, 


Te [GeV] 


102 


102 


103 


103 


T, - Te [GeV] 


0.05 


0.02 


0.02 


0.01 


T, - n [GeV] 


0.63 


0.28 


0.36 


0.05 




0.38 


0.25 


0.25 


0.23 




1396 


409 


624 


187 


F/T 


240 


223 


232 


221 


Ftw/T 


143 


143 


143 


143 


R [GeV^i] 


1.6 


2.9 


2.4 


4.3 


Rli 


8 


10 


9 


6 


d/R 


0.30 


0.25 


0.28 


0.29 



Table 1: Observables of the first-order phase transition for four different 
effective potentials, a Higgs mass of 70 GeV and a top-quark mass of 140 
GeV. 



the coarse-grained free energy in condensed matter physics. An important 
aspect of this approach is that scalar fluctuations are only computed around 
the stationary points $ = and $ = $ of Fa[$,T] and not, as usually 
done, around unstable homogeneous scalar background fields. Hence, the 
perturbative approach is consistent and does not break down due to infrared 
divergencies or negative scalar mass terms. The decay rate (^) is similar 
to the result obtained for scalar electrodynamics [§. The difference in the 
pre-factor is due to the different global symmetries which are spontaneously 
broken. 

Let us finally consider the cosmological phase transition. A rough esti- 
mate of the temperature Tg at which the phase transition ends, is obtained 
by requiring 

TiQtt ~ 1, (90) 



28 



where t ~ 0.03 ■ rripi/T'^ . As an example we choose for Higgs boson and 
top quark masses the values mn = 70 GeV and rrit = 140 GeV. The W- 
boson mass is myy = 80.6 GeV. These masses correspond to the coupling 
constants A = 0.04, ft = 0.80 and g = 0.66. From our discussion in sect. 5 
we know that for these parameters higher order corrections are under control 
so that we can still say that the phase transition is first-order. From the 
free energy F (cf. eq.(^)) we can compute the critical temperature T^, the 
barrier temperature T^, the surface tension a, the correlation length of the 
symmetric phase ^ = and the value of the Higgs field inside the droplet, 
i.e. ipc/Tc. From eqs. (^) and ( pO] ) we then obtain the temperature Tg 
and the size of R of the critical droplet. All these quantities are listed in 
table 1. Note, that the size of the pre-factor in the total rate relative to the 
exponential is of order 1%. The thin wall approximation for the Higgs mass 
used in the table is marginally applicable. We have plotted V on fig. ITl for 
three different Higgs masses at the corresponding temperature Tg. As it can 
be seen, the larger the Higgs mass the better the thin wall approximation. 

How reliable are these results? Since the ratio ^Pc/Tc is rather small, the 
longitudinal plasma mass mi (cf. eq. (p9D ) is essentially independent of (p. 
Hence, the longitudinal degree of freedom of the vector field decouples and 
does not contribute to the tunneling process. This leads to a reduction of 
the cubic term in the effective potential by | 0. This effect is not included 
in the potential V. Let us denote the potential with the reduced cubic term 
by V2. As shown in table 1, the effect on the surface tension is considerable. 
However, the phase transition clearly remains first-order with some change in 
the relevant temperatures and droplet properties. In order to illustrate the 
possible effect of higher order corrections we have also computed the observ- 
ables for the potential including the scalar loops (Va) and for the potential 
with scalar loops and a non- vanishing magnetic mass (7 = 1) (cf. eq. (pTD). 
Note in particular the effect of the magnetic mass on the surface tension. A 
more detailed quantitative discussion will be presented in a separate paper 

To conclude, we have obtained a consistent description of a cosmological 
first-order electroweak phase transition for values of the Higgs boson and top- 
quark masses which are compatible with present experimental limits. The 
phase transition is only weakly first order. The quantitative description of 
the transition can be used as input for models of baryogenesis. 
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7 Summary 



In the previous sections we have studied the transition from the symmetric to 
the broken phase in the SU{2) gauge theory at finite temperature. We have 
seen that, due to infrared divergencies, ordinary perturbation theory to any 
finite number of loops does not yield a useful approximation to the effective 
potential. However, an improved perturbation theory, which takes plasma 
masses into account, describes consistently the symmetric phase {cp — 0) and 
also the broken phase ((/? > 0) in the neighbourhood of the second non-trivial, 
local minimum of the effective potential. Using this improved perturbation 
theory we have evaluated the effective potential including all terms cubic in 
the couplings and shown that all contributions linear in (p cancel. To this 
order in the couplings the ring summation and the improved perturbative 
approach are equivalent. 

The plasma masses have been determined from a set of one-loop gap equa- 
tions. A non-vanishing transverse gauge boson plasma mass was obtained, 
originating from the non-abelian gauge interactions. Based on the gap equa- 
tions we also found a range in the couplings g, X and ft, the temperature 
T and the scalar field (p, where the perturbative approach is reliable. The 
dependence on the top-quark Yukawa coupling ft turns out to be irrelevant. 
Knowing this range in T and (p as function of g and A where the effective 
potential is rehable has allowed us to determine the range in A where the 
symmetric phase is metastable. As a criterion we required that at the origin, 
(p = 0, the effective potential has only a local and not a global minimum for 
the allowed values of T. 

We find that the electroweak phase transition is weakly first-order and 
that our perturbative approach in fact breaks down for Higgs masses close 
to the present experimental lower bound. Higher order corrections from 
scalar loops become important for Higgs masses around 80 GeV, for vanishing 
magnetic mass of the W^-boson (cf. sects. 4A and 4C). For non- vanishing 
magnetic mass (7 = 1) the perturbative approach breaks down around mn = 
70 GeV. 

Following the theory of Langer we have finally computed the nucleation 
rate for critical droplets, and we have discussed some aspects of the cosmo- 
logical phase transition. Up to Higgs masses of order 80 GeV the picture 
of a first-order transition, which proceeds via nucleation and growth of crit- 
ical droplets, appears self-consistent, if the magnetic plasma mass vanishes. 



30 



However, the thin wall approximation is only marginally applicable. Fur- 
thermore, the Higgs vacuum expectation value inside the critical droplet is 
much smaller than the value required by models of electroweak baryogenesis. 

Our results could be improved in several respects. Clearly, a complete 
computation of the effective potential to order g'^, and would be very 
valuable in order to further test the convergence of the perturbation theory. 
A method to identify the relevant contributions is given in appendix A. Fur- 
thermore, the validity of the expansion in powers of derivatives used in sect. 
P has to be examined in greater detail. However, the most crucial ingredient 
concerning the order of the transition and the difference between abelian and 
non-abelian gauge theories appears to be the magnetic gauge boson mass, 
whose origin and size require further investigations. 

We would like to thank M. Liischer, I. Montvay, M. Renter, N. Tetradis 
and C. Wetterich for helpful discussions and comments. 
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Appendix A Ring diagram summation 



As it was promised we will now prove that only ring diagrams give non- 
zero contributions cubic in the couplings {g, a/A and ft). 

Instead of g, ^/X and ft we will use in the following a generic coupling h. 
The naive order of a Feynman diagram is given by 

naive = ^naive ~^ '^^naivei (■^•l) 

where V^J^^^ is the number of the vertices of order h (the triple vector ver- 
tex, the ghost-ghost-vector vertex, the scalar-scalar-vector vertex and the 
fermion-fermion-boson vertices) and V^^^^^ is the number of the vertices of 
order (the scalar triple self-coupling, the vector-vector-scalar vertex and 
the quartic boson vertices). We study a graph with non-zero number of ver- 
tices contributing to the effective potential. It is one-particle irreducible with 
no external lines. Momentum conservation gives a relationship between the 
total number of vertices {V = V^a]?;e + KlaiLe); the number of independent 
loops (L) and internal lines (/), 

V + L-I = l. (A.2) 

The momenta of the internal lines qi,q2, ■■■,qi are linear combinations 
of the loop variables pi,p2, ■■■,Pl- A lower order contribution in h {O < 
Onaive) could arisc if some of the Matsubara frequencies, e.g. Q'lo, '?20; •••Q'sO) 
vanish. Denote the loop momenta with vanishing Matsubara frequencies 
as pi,p2, --jPr- We call lines or loops with vanishing Matsubara frequencies 
soft, those with no n- vanishing Matsubara frequencies hard (cf. ref. ||30|| ). 
Fermions have always non-zero Matsubara frequencies, thus fermion lines 
and loops are always hard. 

The substitution in the soft loop variables 

Pi = mi^^yi , p2 = mf''y2 , ... , Pr = ml^^yr (A. 3) 
results in a similar substitution for the propagators 

qi = rn'pxi , q2 = 1112^X2 , ... , (fs = m^Pxs. (A.4) 

Here the masses are proportional to the coupling h 

if/ 2 for vector boson; 
mf = haf' or I) ^ ai= \ ^ifi'^ — v'^ for Goldstone bosons; (A. 5) 

V3(p^ — f 2 for the Higgs boson . 
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The above transformation gives for each soft loop integration, 

J d^pi = ml j d^Vi (xh^ J d^Vi , (A.6) 
an additional factor of h^, and for each soft boson propagator, 

= -^D{xi) oc h-^Dixi) , (A.7) 

an additional factor of h^"^. The same factorization can be done for vertices 
with momentum dependence if all of the incoming lines are soft, 

y\^{Pi,Pj,Pk) oc hyV{xi,Xj,Xk) , (A. 8) 

where W is a linear combination of the soft momenta. Note, that for the 
ghost-ghost-vector vertex this additional factor h appears not only in the 
case where all of the lines are soft, but also if one of the incoming ghost lines 
are soft. 

Hence, inspection of the Feynman rules suggest that the diagram is at 
least of the order 

O = 2^(2) + V^^^ + 3r - 2s, (A.9) 

where V^^^ is the number of the vertices with naive order h and at least two 
incoming hard lines. 

Remove now the r soft lines from the graph and consider the remaining 
subdiagram containing only hard lines. This subdiagram is not necessarily 
connected but all of the j = 1, ...,n connected parts are closed graphs with 
no external lines. (The original graph has no external lines and due to the 
momentum conservation at the vertices there is no way to connect a hard line 
only to soft lines.) The relationship ([A.2|) is valid for each of these connected 
graphs with Vj vertices, Lj loops and Ij lines, thus 

V, + L,-/, = 1, j = l,..,n . (A.IO) 

Note, that in this case some of the vertices are connected only to two lines 
but this fact has no influence on our consideration. Clearly, the number of 
the vertices in the j-th connected graph Vj is the sum of the number of the 
vertices of order h and the number of the vertices of h"^, thus 

Vj = vl'^ + Vl^l (A. 11) 
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The sum of the hard hnes in the individual connected hard sub diagrams is 
the total number of the hard lines (/ — s), and the sum of the individual hard 
loops is the total number of the hard loops (L — r), 

j j 

Since the vertices of order h have at least two hard lines the connected hard 
subdiagrams contain all of them 



^V,.W =01). (A. 13) 

j 



Summing ( [A.10|) over j one obtains 

J2{V,-l) + {L-r)-{I-s) = 0. (A.14) 
j 

Solving for s and inserting it into ( |A.9| ) yields 

O = 2V^^^ + 2V^^^ - V^^^ +r + 2r-2s 

= 2V - V^^^ +r + 2r + 2\ L-r-I- Y,iVj - 1) | • (A.15) 



From this equation, the relation V + L — I = 1, ( A. 11 ) and (|A.13 ) one obtains 
the final answer 

= {2 + r) + + ^ - 2). (A. 16) 



The case r = corresponds to no soft lines, thus the order of the diagram 
is just the naive order. For r = 1 the order of the diagram is at least 3, 
where we have one soft loop. If V^''^'' + 2Vj'^'^ — 2 = for j = 1, then 
the order of the graph is still 3. To ensure this for a given j there are two 
possibilities. Either there is only one quartic coupling of order h"^ in this hard 
subgraph or there are two triple vertices of order h (since we only study one- 
particle irreducible contributions to the effective potential it is not possible 
to have only one triple vertex in a hard subgraph.) These are just the ring 
diagrams. Note, that the self-energy insertion with two scalar triple vertices 
or vector-vector-scalar vertices do not belong to this class. 
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If one wants to collect the graphs of order h'^ one has to take into ac- 
count the previous (r = 1) graphs one order further in the high temperature 
expansion, a finite number of graphs with no soft loops (r = 0), no new 
terms with one soft loop {r — 1) and the above mentioned hard graphs 
+ 2V^'^^ -2 = 0) for the case with two soft loops (r = 2). As before, 
self-energy insertions to order are satisfactory. 
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Appendix B Partial summation of superdaisy graphs 



As we have shown performing an iterative solution of the gap equations 
one can obtain higher order contributions to the effective potential. In this 
appendix we calculate the g^, and contributions. Clearly, not all of 
this corrections are involved in the solution of the gap equations, because a 
summation where all the self-energy insertions are calculated at zero external 
momenta cannot be justified. 

The leading order masses ( |5^ were obtained by one iteration of the gap 
equations, keeping the terms of order g, a/A or ft. Inserting these masses 
in the gap equations and performing a Taylor expansion up to order g"^, A 
and ft in the square-roots of the obtained quantities, one gets the next order 
masses. 

= ga+(-^-^-^ ^]g^T.-\ 

^ ^ '32 a 32 a 2 a 32a(c + a)r 



(2) 




dc 



96 6 (c/ + c) 48 6 (c/ + c) 

2t^„-1 



c2 



1 Tti^ 
96 6 (rf + c) 48 6 (c + 6) ; ^ 

3y?2 3 5 3 a 3^,2 




64 c6 16 c 32 c 128 ca 

3 (i , \ 1 

3/8 TX-R-^ 

8c ' j 

A i TT g , 



%dc 8 c2 



2,^2 



TAV' 



2d^% (rf + c) 



(B.l) 



In order to reduce the size of our formulas a, 6, c and d were introduced. 
They are proportional to the leading order masses. Thus, 



(1) 
mV 
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b = 



d 



m 



(1) 



m 



(1) 



m 



(1) 



9 



2 



A 



A 



yi6 2^2 4^ 



1^ 



(B.2) 



By use of this masses the next iteration with a Taylor expansion of order 
g'^, A^ and can be carried out. After this straightforward but fairly tedious 
calculation the next order Goldstone mass term is of the form: 



(3)2 / 'P 3V?^ \ 4 2 -2 -4 

tTT'i = : :^ : 9 : 9 — A i tt q 

^ \ 2{d + cfd 8{d + cf c2 8 (d + c)^ cd; 



' 8(d + c)'d (rf + c)' 8ci((i + c) 



9^ _ 3^2^; -3^2^-2-2 

+ 32c2 S{d + cfc ?>2cd)^' ^ 

V64d 128 c 64 c 256 c6 512 ca 128 rfy ^ 
rA2^^ ^ {-l-'-i)\Tg ^ (-!f-f)V 



+ 



g^TT ((i + c) TT TT 

^ rf2 36 9rf c2 



^ 256 6(ci + c) 32 a 512 a l2Sh{d + c) 256 6(d + c) 
3 c 3(^2 ^ 



, _| — r I r I ji2 4 -2 

512 a 512a(c + a) 128 6(c + 6) j 



32 (d + c)c 16(d + c)^d 128 (rf + c)ca 



32{d + cfd 6A{d + cfcb 16 {d+cfc^ 32 d 
3d 



+^lr^AV-. (B.3) 
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Similarly, it is easy to get the corresponding Higgs mass-squared of order 
and 

^ 512c3a+ / 64^3 128 rf3 

27(^2^ 27v?2^ 9rf 3c 27y?^ \ 
~ 128 c3 ~ 64 c3 ^ 3215 ^ 32rf ~ 256 c^fo j 

+ — ^ H ^ : H ^ \ : A T^g^-K'^ 

V256c6 128 c 64 d 64 c 512 ca 128 

V 4c 4rf y 

, / ^ _ 27^'^ 15^\ ,3^2 -2 -2 

1^ 32rf3 32c3 8rf(d + c) 32cd ?,2(Pj ^ 

( 3 b £ cd 

+ . . . ^ + tt;^ - ....... ^ + 



256 6(rf + c) 32 a 256 6 (c/+c) 128 6((i + c) 
^ ^ 9 c/ 3y,25 3(^4 9^2^ 



128 6(c + 6) 512 a 128 a^ 2048 a3(c + a) 2048 a^ 
3v?^c (p^d^ ({f^ (p^dc 



+ 



2048 a3 1024 63(c/ + c) 512 63 (c + 6) 512 63(c;+c) 



3^' ^'c' 3 c 



2„4_-2 



512a(c + a) 1024 63(c; + c) 512 a 
V 64 a 32 6 8 16 / ^ 



81 (/?^ 3(/?^ 21 



_|_ I " r -■ T 1 ■v4Tn2 -2 -4 

1^ 32 c4 8rf3(d + c) 32c3rfy' ^ 

TT 

According to (|49|) the masses determine the effective potential of order g^, 
A2 and 

1^(3)(<^,T)= r d^'mf\^\T). (B.5) 



X 

Due to the complicated structure in ip of the integral, (a, c and d are irrational 
functions of if) we have not evaluated it analytically. Instead of it, we have 
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used a numerical integration in order to determine the effective potential at 
the critical temperature and to calculate the surface tension. 
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Appendix C Coset space volume 



Due to the global SU (2) x U (1) symmetry of the standard model, which is 
spontaneously broken to the electromagnetic U{1) subgroup, the fluctuations 
around the saddle point (cf. eq. (|79|) ) contain three zero modes which have 
to be treated in the usual way by the method of collective coordinates. 

An inflnitesimal transformation of the saddle point $, which yields an- 
other saddle point with the same free energy, is given by 

3 

= dU^ = i duJiTi^ , (C.l) 

i=l 

where 

7^1,2 = ^n,2, T^ = ]^{t^-E) ■ (C.2) 

here E denotes the identity matrix. Transformations generated by T4 = 
^(ts + E) leave the saddle point $ invariant. Comparing eq. ( |(J.1D with 

a \/-\ -\- Id \/!~, 1 

(C.3) 



yields the connection between the flelds dxi and the group parameters dui. 
The corresponding measure of the functional integral is 

[dx] = n (^) ' dc,^ (C.4) 

where dxi = ddci, and is a normalized function, i.e., / d^x^f = 1. Eqs. 
3), O and O yield 




[dx] = --^l^-^ldWyt[du;.. (C.5) 
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From eqs. (|C.1|) and (|C.2|) one can easily calculate the metric on the 
group manifold at the origin a; = 0: 

g,, = -tr{U-'—UU-'—U) = c(z)5,„ (C.6) 



40 



where 

c(l) = c(2) = i c(3)=c(4) = l, (C.7) 



detg{u = 0) = - . (C.8) 

Globally, the symmetry group acting on the Higgs field is not SU{2) x f/(l), 
but rather 5*0 (3) x Z2 x 5*0(2), where the elements of Z2 are E and —E, and 
5*0(2) corresponds to phase transformations with the phase varying from 



to TT. For the normalization (|C.8|) the volume of this group is (cf. ref. |31]) 

J duji...dujiy/g = 27r^. (C.9) 

The Higgs vacuum expectation value breaks this symmetry to the subgroup 
[/(I). The corresponding volume of the coset space is 

V{S0{3) X Z2X S0{2)/U{1)) = vrl (C.IO) 

For the integration over the zero modes one now obtains 

The corresponding three eigenvalues in the symmetric phase are ~ /i^, where 
/i = m^(0,T). Hence, the relative contribution of the zero modes to the 
nucleation rate F is //^ / [dx] ■ 
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Figure caption 

Figure |l] : One-loop contributions to the effective potential. 

Figure ^ : Two-loop contributions to the effective potential. 

Figure ^ : Ring diagram contributions to the effective potential with one- 
loop self-energy insertions. 

Figure § : One-loop self-energy contributions to the gauge-boson. 

Figure ^ One-loop contributions to the effective potential with full propa- 
gators including counterterms. 

Figure ^ The gap equations: All one-loop self-energy corrections with full 
propagators. 

Figure ^ The maximal value of the Higgs mass for which the phase tran- 
sition is first-order as function of 7. 

Figure ^ The maximal Higgs mass as function of the convergence param- 
eter ^. (The dependence on the top-quark mass is smaller than the 
width of the line.) 

Figure ^ The effective potential to order and A'^/^ at the smallest tem- 
perature (T*) allowed for ^ = 2. 



Figure |T^: Surface tension as function of the Higgs mass for different values 
of 7. 

Figure O: Surface tension as function of the Higgs mass with the 

contributions (full line), with A-contributions only (long-dashed line) 
and with contributions from the partial summation of and A^ terms 
(short-dashed line). 

Figure |l^: The potentials U^{r) and U^{r) for the scalar fluctuations (p and 
X- 

Figure |1^: The potential V for three different Higgs masses at the corre- 
sponding nucleation temperatures Tg. The potential for mn = 100 GeV 
(200 GeV) has been multiplied reduced by a factor 10 (factor 400). 
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